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Abstract 

We study the curvature strength and visibility of the central singularity arising in Tolman-Bondi- 
de Sitter collapse. We find that the singularity is visible and Tipler strong along an infinite number 
of timelike geodesies, independently of the initial data, and thus stable against perturbations of the 
latter. 
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1 Introduction 

It has recently been shown by Goncalves |EJ and, independently, by Deshingkar et al. that a central 
curvature singularity develops in asymptotically de Sitter spherical dust collapse. When curvature 
growth is examined along radial null directions, the singularity is found to be visible, provided the 
initial central density distribution, p c {r), obeys certain differentiability criteria p]]. Let n € N + be 
such that (d n p c (r) / dr n ) r= o = p n is the first non-vanishing derivative of the central energy density 
distribution. The singularity is visible for n < 3. For the cases n < 2, a single radial null geodesic is 
emitted from the singularity; for n = 3, if, in addition, p% obeys certain constraints, an entire family of 
outgoing radial null geodesies has its past endpoint at the singularity. Regarding curvature strength, the 
singularity is Tipler strong Q| for n = 3, and may or may not be Tipler strong for n < 2. This behavior 
is analogous to the one found in asymptotically flat spherical dust collapse, whose analysis along radial 
null directions shows the central singularity to be visible for n < 3 j|, ||, |(|, and gravitationally strong 
for n = 3 (weak for n < 2) |J. 

All the previous results were derived based on an analysis of curvature growth along radial null 
geodesies. However, such analyses, based solely on radial null geodesies, are too restrictive and have 
given rise to the misleading notion that the physical nature of the singularity in asymptotically flat dust 
collapse is critically determined by the initial data. A notable exception is the work of Deshingkar, Joshi 
and Dwivedi 0; as their analysis shows, when one considers timelike radial geodesies, the singularity 
is found to be locally naked and Tipler strong for an infinite number of non-spacelike radial geodesies, 
irrespective of the initial data. A recent detailed analysis of non-spacelike geodesies in (asymptotically 
flat) dust collapse by Deshingkar and Joshi || showed that a family of radial null and timelike outgoing 
geodesies has its past endpoint at the singularity, thereby confirming the results of |Q . 

In this paper, we perform an analysis of visibility and curvature growth along radial timelike direc- 
tions for the case of Tolman-Bondi-de Sitter (TBdS) collapse. We also find that the central curvature 
singularity is at least locally naked, and Tipler strong, independently of the initial data, which implies 
stability against perturbations of the latter. Such a behavior is not unexpected, since the structure of 
the singularity in TBdS collapse is quite similar to that of asymptotically flat dust collapse [j], 0, ^, 0, . 
This result adds further evidence to the robustness of dust collapse against the introduction of a cos- 
mological constant. 

Geometrized units, in which G = c = 1, are used throughout. 
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2 Central naked singularity in Tolman-Bondi-de Sitter collapse 

The Tolman-Bondi-de Sitter family of solutions is given by a spherically symmetric metric, written here 
in comoving coordinates {t, r, 9, </>}: 

ds 2 = -dt 2 + -dr 2 + R 2 dQ 2 , (1) 

1 — k 

where ' = d r , together with the stress-energy tensor, 

T ab = p(t,r)S t a S t b -^-g ab , (2) 

where p(t, r) is the energy density distribution, and A > the cosmological constant. In spherical 
symmetry there are no gravitational degrees of freedom, and the metric is thus uniquely determined 
from the initial data, p(0, r) and R(0, r), where ' = d t and R(t, r) is a solution of 

R 2 = — -k + -R 2 , (3) 

with initial data given by 

m'(r) = 4TrR 2 R'p{t,r), (4) 

where m(r) is the Misner-Sharp mass p]| (which agrees with Hawking's quasi-local mass Jl2[ and 
ADM mass, in the appropriate limits) and is uniquely determined from the initial density profile, p(0, r). 

The case k = (corresponding to gravitationally unbound matter configurations) is qualitatively 
equivalent to those with k ^ 0, and it allows for an analytical solution of Eq. (^) in closed form: 

R(t,r) = (^ 3 sinh 2 / 3 T(t,r), (5) 

T(t,r) = ^[t c (r)-t], (6) 

where t c (r) is the proper time for complete collapse of a shell with initial area radius R(0, r), which is 
fixed by Eq. (|) at t = 0: 

where the scaling i?(0, r) — r was adopted. 

The relevant derivatives of the area radius are 

R'(t,r) = H^ + yi^cothTj, (8) 
R(t,r) = -^RcothT. (9) 
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At t = t c , R(t c ,0) — 0, and the Kretschmann scalar 

l j ml mm! m' 2 8 , , , , . 

JC = i? Qbc(i i? afctrf - - 8^ + 12-^ + -A(A + 4^), (10) 



diverges, thereby signaling the existence of a curvature singularity at r — 0. This singularity will be at 
least locally naked if the outgoing radial null geodesies equation, 

^ = R>, (11) 
ar 

admits a regular solution, for given initial data, at t = t c (0) — to. Expanding p(0, r) = p c {r) near r = 0, 

PcM=]>>r l , (12) 

near the singularity (where r, T — > + ) we have, to leading order, 

m(r) = m r 3 + m„r™+ 3 +0(r n+i ), (13) 
tc(r) = < + ^"+O(r" +1 ), (14) 

1 

R(t,r) = (mJr + M„r n+1 )(i + t„r n -t)$ + 0(r"+ 2 ) x O(tI), (15) 

where t n and M„ are real coefficients linear in m n = (47r/n)p n , with n > 0; p n = {d n p c /dr n ) r= o is the 
first non-vanishing derivative of the central energy density distribution, and 



t = sl nh-(yA.j. (16) 

Assuming that there is a regular solution to the outgoing null geodesies equation near r = 0, which 
we write, to leading order, as 

i = i +ar CT , (17) 

where a, a 6 M + , the singularity will be at least locally naked if Eq. flll] ) admits a self-consistent 
solution, and covered otherwise. Since the geodesic must lie on the spacetime, from Eqs. (|l4|) and ( |l7[ ) 
it follows that a > n. From Eqs. (Uq) and (E7|) , we obtain 



R(t,r)= f^Vt|r* +1 +0(r CT + 2 -t). (18) 



V 2 

Let us first consider the a > n case. Differentiating Eqs. ( p7| ) and ( |l8| ) w.r.t. r, and requiring a 
self-consistent solution of Eq. ( pd] ) yields 

<7 = 1 + (19) 
a = {^ftl (20) 



The condition a > n now reads n < 3. For n — 1, 2 (i.e., for pi 7^ 0, or p\ — and pi ^ 0) there is a 
self-consistent solution to the outgoing radial null geodesies equation in the limit t — > to, r — > 0, and 
thus there is at least one outgoing radial null geodesic starting from the singularity. 
For n = a = 3, an outgoing radial null geodesic exists, provided [|l] 

— - 4 + VM 4 + 16 - AM 2 < i 3 < -M 4 (21) 
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where M = (9mo/2)3. For the special case (with n = 3) when M > 8.799799016 and ^3 is uniquely 

determined from po [cf. Eq. (45) in M], an entire one-parameter family of outgoing radial null geodesies 

departs from the singularity. Clearly, this is a set of measure zero in the initial data. 



3 Visibility 

Let us consider radial timelike geodesies (RTGs), with tangent vector K a — where t is an affine 
parameter along the geodesic, and 



K l = ±^J\ + R' 2 (K r ) 2 , (22) 

K r (K r ) 2 
K r R' + 2K r R' + — (K r )'R' + ±-4-R" = 0, (23) 

where the first equation is simply K a K a = — 1, and the second follows from the geodesic equation. By 
inspection, one sees that Eq. ( p3| ) admits the trivial solution 

K l = ±1, (24) 

K r = 0, (25) 

which leads to 

t = t ±(T-T ), (26) 

r = ro = const., (27) 

where the plus or minus sign refers to outgoing or ingoing RTGs, respectively. The outgoing RTG 
departing from the singularity is given by r — 0, and t = to + t — tq, and thus does not belong in the 
spacetime. The ingoing RTG is given by r — 0, t = to — r + To, where to = t c (0) = is the time at 
which the RTG arrives at the singularity. 

3.1 Genericity 



Since Eq. (|23|) is a mixed first-order linear PDE for K r (t,r), its solution need not be unique |13j. 
Indeed, one can explicitly show that there are infinitely many other solutions, as follows. Near the 



singularity, we can write, to leading order 

*rtgM = t Q + br p , (28) 

R' = a ir \ (29) 

R' = a 2 r q - p , (30) 

R" = a 3 ri-\ (31) 

where b,p,q G R + , and a, € K (i = 1,2,3) are constants uniquely determined from the initial data, 
and the t dependence of R has been absorbed in the exponent q [which can be done without loss of 
generality, since R is evaluated along the geodesic defined by Eq. (p8|)]. 

Let us now assume that K r (t, r) oc (t - t ) a r f3 , where a, j3 G K. From Eq. (f|) 

wc have then 

#r = Cir ap+/5 ; (32) 

K r = c 2 r {a - 1)p+(} , (33) 

(JT)' = csr"^- 1 , (34) 



where c, G M (i = 1, 2, 3), are free constants, yet to be determined. With these ansatze, Eqs. (|22|)-([23|) 
read: 



K l = ±y/l + a 2 c 2r 2 ('?+«P+/ 5 ), (35) 

aiC2r (^i)p+/3-9 + 2a 2Cl r( Q - 1 ^-« + ^™ r 2(a P +/3)+ 9 -i + ^ r 2(a P +«+ g -l = Q (36) 

If the free constants Cj and p,q,a,/3 are such that the system (p2|)-(p3|) is self-consistent, then there 
are outgoing radial timelike geodesies emanating from the singularity — each one of them uniquely 

determined by a particular set {p, q, a, f3,Ci,i — 1, 2, 3}. We now look for solutions of the form K l — 
const., and in which the left- hand-side of Eq. ( |3^ ) is homogeneous in r. This implies 

ap + q + (3 = 0, (37) 

ap + q + (3-p = 2(ap + (3) - 1 + q, (38) 

which is solved by 

1-/3 

P=l + q=—^. (39) 
1 + a 

Equation ( p3| ) becomes then 

C(c 4 )rfc^ = 0, (40) 
which has an infinite number of solutions, given by the algebraic constraint 

C(ci) = aic 2 + 2a 2 ci ± -===- = 0. (41) 



Hence, there are infinitely many radial outgoing (ingoing) timelike geodesies, with past (future) endpoint 
at the singularity. The existence of these families of RTGs is independent of the initial data — which 
enters the algebraic constraint via the at constant coefficients — since the algebraic constraint C(cj) = 
always has a solution, due to the three degrees of freedom for the Ci coefficients. 

One can also prove the existence of an infinite number of timelike geodesies with endpoint at the 
singularity, by resorting to simple causal structure considerations, as follows. 

The existence of a future-directed timelike geodesic cr(A) with future endpoint at the singularity, 
implies that its chronological past I~ [cr] is a timelike indecomposable past set (TIP) p] . In other words, 
the TIP is generated by cr, which is timelike and future-inextendible. The locally naked singularity 
itself — which is the future endpoint of cr, er(Ao) — constitutes a singular TIP which contains the past 
I~(p) of any point p £ ct\{(t(tq)} fjsfl . Consider now another point q a in the past of the singularity 
and in the chronological past of p. Since I~ (q) C I~ (p) , it follows that there exists a timclikc curve from 
p to q, say £, satisfying C Consider then a small compact ball B — with a suitably defined 

"radius", e.g., proper geodesic distance along cr, d CT (A, Ao) — in the neighborhood of o~(tq), partially 
contained in TIP I~ [cr] and in the chronological future of p, i.e. , B n I~ [cr] = C and I + [C] C I + (p) . 
For each point x G C, there is a timelike geodesic from p to x. Since C is compact, there is an infinite 
number of such points that can be joined by timelike geodesies from p € a. Now, let p be an arbitrary 
point on cr and take the limit of approach to the singularity, do- (A, Ao) — * + ; it then follows that there 
is an infinite number of future-inextendible timelike geodesies with future endpoint at the singularity. 

3.2 Apparent horizon 

The existence of outgoing timelike geodesies with past endpoint at the singularity renders the latter 
visible if and only if the geodesies are emitted before or at the time of formation of the apparent horizon 
(AH). The AH is a spacelike 2-surface defined by the outer boundary of a trapped region, and in the 
adopted coordinates it is given by R, a R^ b g ab = 0. With the metric ^ and Eq. (g) we have then 

-R 3 - R + 2m = 0. (42) 
3 

This equation has three distinct real roots if 3my/X < 1, two of which are positive and given by 

(43) 



2 

R± = —= sin 

VA 



- sin _1 (3mVA) + ^-5(1 ± 1) 



with i?_ > R + > 0, corresponding to the choice < sin u> < tt/2, < uj < 1. The third root, R, 
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— R_ — is negative and hence unphysical. R- is a generalized cosmological horizon (i?_ = ^/3/A, 
when m = 0) and R + the black hole apparent horizon (i? + = 2m when A = 0; the apparent and event 



horizons coincide in the static case). For 3mvA = 1, the two horizons coincide. If 3mvA > 1, there is 
one negative real root and two complex (conjugate) roots, all of which are unphysical: the spacetime 
does not admit any horizons in this case. From the "+" solution in Eq. (f43|), together with Eqs. (|])-(@), 
we obtain 

9 . f / S r /I . ~- \ 1 3/2 1 

(44) 



2 

tAa(r) = t c (r) == sinh" 1 



sin f — sin (3mvA) 



VA \ V 3mA 

At the origin we have t AH (0) = t c (0)-£, where < £ < (2/%/3A) sinh^ 1 ^V3/(mA). Hence, £ O rtg(0) > 
^ah(O). If ioRTc(O) > £ah(0), any radial geodesic emitted from the singularity at r = will do so after 
the AH forms, and is therefore unavoidably trapped — the singularity is covered. If ^ortg(O) = £ah(0), 
the singularity is at least locally naked, provided (^ a ) r=0 > ( dt °£ TG ) r=0 - For £ = 0, ^ml = ^ > 
(i.e., there are no shell-crossing singularities); if p ^ 1 in Eq. (|2S|), then ( dto J^ TG ) _ = 0' an d the 
singularity is thus always visible. For the special case p = 1, appropriate choice for 6 in Eq 
renders the singularity visible. 



4 Curvature strength 

Let 7(7") be a RTG, with tangent vector K a — that terminates at the singularity at r = To. Such 
a singularity is said to be Tipler strong along 7(7") if the volume three-form V(t) vanishes in the limit 
T ~~ * r o |§- If the scalar = R a bK a K b obeys the strong limiting focusing condition, 

lim (r - T ) 2 * > 0, (45) 

then the singularity is gravitationally strong in the sense of Tipler \V^ . This sufficient condition 
guarantees that any three-form defined along 7(r) vanishes in the limit r — > to, due to unbounded 
curvature growth. 

From Eqs. (0)-(|) and (§|)-©, we have, 

* = A A. (46) 

R 2 R> y ' 

Using Eqs. and ©, we obtain, near r = 0, 

* = £»T* (V + ^*«r n ?7-^ , (47) 
T) = t n r n T(r-r Q ). (48) 



Hence, 



lim(r-r„) 2 f= |, (49) 

T— >T<] O 
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and the singularity is therefore Tipler strong. This is a rather robust result, in that it is independent of 
the initial data and it holds true along an infinite number of timelike geodesies. We note that, since all 
the sufficient criteria for Tipler strong singularities are sufficient for deformationally strong ones Jt8| , |l9| , 
the central curvature singularity in TBdS collapse is also deformationally strong. 



5 Concluding remarks 

We have shown that the central singularity arising in spherical dust collapse with a positive cosmological 
constant is Tipler strong and at least locally naked, along an infinite number of timclike geodesies, 
independently of the initial data. Such a singularity could possibly constitute a counter-example to 
the strong cosmic censorship conjecture [fl5| , [HI , if it were to remain stable against perturbations of 
matter and geometry. The work of Harada et al. J2(| and Iguchi et al. [2^], showed that the central 
curvature singularity in asymptotically flat dust collapse is marginally stable against non-spherical linear 
perturbations. These authors examined the dynamical evolution of non-spherical metric, matter, and 
matter coupled to matter perturbations in Lemaitre-Tolman-Bondi collapse. They found that whilst the 
singularity itself is not a strong source of gravitational radiation, the non-spherical perturbations remain 
well-behaved and finite in the vicinity of the singularity and in the limit of approach to the Cauchy 
horizon. Given the similar structure of the central singularity in Lemaitre-Tolman-Bondi collapse and 
the one discussed in this paper, one may speculate that similar stability properties also hold when 
A > 0. This issue is currently being investigated. 
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